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Abstract. We investigate the possibility approximating relativistic effects in hydrodynamical simulations of stellar 
core collapse and post-bounce evolution by using a modified gravitational potential in an otherwise standard 
Newtonian hydrodynamic code. Different modifications of a previously introduced effective relativistic potential 
are discussed. Corresponding hydrostatic solutions are compared with solutions of the TOV equations, and hy- 
drodynamic simulations with two different codes are compared with fully relativistic results. One code is applied 
for one- and two-dimensional calculations with a simple equation of state, and employs either the modified ef- 
fective relativistic potential in a Newtonian framework or solves the general relativistic field equations under the 
assumption of the conformal flatness condition (CFC) for the three-metric. The second code allows for full-scale 
supernova runs including a microphysical equation of state and neutrino transport based on the solution of the 
Boltzmann equation and its moments equations. We present prescriptions for the effective relativistic potential 
for self-gravitating fluids to be used in Newtonian codes, which produce excellent agreement with fully rela- 
tivistic solutions in spherical symmetry, leading to significant improvements compared to previously published 
approximations. Moreover, they also approximate qualitatively well relativistic solutions for models with rotation. 
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1. Introduction 

It is a well known fact that gravity plays an important 
role during all stages of a core collapse supernova. Gravity 
is the driving force that at the end of the life of massive 
stars overcomes the pressure forces and causes the collapse 
of the stellar core. Furthermore, the subsequent super- 
nova explosion results from the fact that various processes 
tap the enormous amount of gravitational binding energy 
released during the formation of the proto-neutron star. 
General relativistic effects are important for this process 
and cannot be neglected in quantitative models because 
of the increasing compactness of the proto-neutron star. 
Therefore, it is important to include a proper treatment 
of general relativity or an appropriate approximation in 
the numerical codes that one uses to study core collapse 

supernovae. 

Recently iLiebenddrfer et all l|2005l) performed a com- 
parison of the results obtained with the supernova simula- 
tion codes Vertex and Agile-BoltzTran which both 
solve the Boltzmann tr ansport equation for neutrinos. 
The Vertex code (see iRampp fe Jankal 12002^) is based 



Send offprint requests to: H.-Th. Janka, 
e-mail: thj@mpa-garching.mpg.de 



on the Newtonian hydrodynamics code Prometheus 
l|Frvxell et alJll990|) and utilizes a generalised potential to 
approximate relativistic gravity. The Agile-BoltzTran 
code of the Oak Ridg e- Basel group l|Liebendorfer et all 
l200ll 120021 12004L 12005(1 is a fully relativistic (ID) hydro- 
dynamics code. The comparison showed that both codes 
produce qualitatively very similar results except for some 
small (but growing) quantitative differences occurring in 
the late post-bounce evolution. Inspired by this compar- 
ison we explored im provements of the effective relativis- 
tic potential used bv lRampp fe Jankal (2002) in order to 
achieve an even better a greement than that reported by 
iLiebendorfer et all <)2005f) . To this end we tested different 
variants of approximations to relativistic gravity employ- 
ing two simulation codes. 

On the one hand purely hydrodynamic simula- 
tions were performed wi t h the code CoCoNuT of 
iDimmelmeier et"al1 l|2002al 120051) assuming spherical or 
axial symmetry. This code optionally cither uses a 
Newtonian (or alternatively an effective relativistic) grav- 
itational potential in a Newtonian treatment of the hydro- 
dynamic equations or solves the general relativistic equa- 
tions of fluid dynamics and the relativistic field equations. 
The latter are formulated under the assumption of the 
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conformal flatness condition (CFC) for the three-metric, 
also known as the Isenberg- Wilson-Ma thews approxima- 
tion lllsenberdll973 IWilson et alJll99d) . which is identi- 
cal to solving the exact general relativistic equations in 
spherical symmetry. This allows for a direct comparison 
of different effective relativistic potentials with a fully rel- 
ativistic treatment in both spherically symmetric and ax- 
isymmetric simulations using the same code. 

On the other hand, we used the computa tionally ex- 
pensive Vertex code |Ramnn fc ,Iankall2002|) for spheri- 
cally symmetric supernova simulations including neutrino 
transport and a microphysical equation of state (EoS). 
The results of these calculations performed with different 
effective relativistic potentials are compared with those 
obtained with the general relativistic Agile-BoltzTran 
code of the Oak Ridge-Basel group (cf. lLiebendorfer et alJ 
120051) . 

In this paper we discuss the results of our test calcu- 
lations and present improved effective relativistic poten- 
tials which are simple to implement without modifying 
the Newtonian equations of hydrodynamics in an exist- 
ing code, and which approximate the results of general 
relativistic simulations very well. 

The paper is organised as follows: In Sect.FJlwe describe 
the effective relativistic potentials used in this investiga- 
tion. In Sect. OH we briefly discuss the most important fea- 
tures of the numerical codes we used for our hydrodynamic 
simulations. In Sect. 0] we present the results obtained by 
applying the effective relativistic potentials in spherically 
symmetric simulations of supernova core collapse and to 
neutron star models, while Sect. is devoted to a dis- 
cussion of the multi-dimensional simulations of rotational 
supernova core collapse. Finally, in Sect. we summarise 
our findings and draw some conclusions. 

Throughout the article, we use geometrised units with 
c = G= 1. 

2. Effective relativistic potential 

Approximating the effects of general relativistic gravity in 
a Newtonian hydrodynamics code may be attempted by 
using an effective relativistic gravitational potential <I> c ff 
which mimics the deeper gravitational well of the rela- 
tivistic case. In the following Sects . |2~D and l2~2l several of 
these effective relativistic potentials will be discussed. 

2.1. TOV potential for a self-gravitating fluid 

For a self-gravitating fluid it is desirable that an 
effective relativistic potential reproduces the solution 
of hydrostatic equilibrium according to the Tolman- 
Oppenheimer-Volkoff (TOV) equation. With this require- 
ment in mi nd and comparing the relativistic eq uation of 
motion (cf. Ivan Riperlll979l iBaron et aljfl9 89) with its 
Newtonian analogon, iRampp fc Jankal l|2002|) rearranged 
the rela tivistic terms into an effectiv e relativistic poten- 
tial fsee lKippenhahn fc WeigerHll99fi for the hydrostatic, 
neutrino-less case). 



Thus for spherically symmetric simulations using a 
Newtonian hydrodynamics code the idea is to replace the 
Newtonian gravitational potential 

$(r) = -47r/ d/r' 2 , 9 ,, (I) 
Jo \r-r'\ 

by the TOV potential 

*rov(r) = -4* jT % (=p + r «(P + p „) 

to obtain the effective relativistic potential $ e ff as 

$cff = $tov- (3) 

Here p is the rest-mass density, e = pe is the internal 
energy density with e being the specific internal energy, 
and P is the gas pressure. The TOV mass is given by 

m T ov(r) = 4:71 J dr' r' 2 (p + e + E + , (4) 

where p„, E, and F are the neutrino pressure, the neu- 
trino energy density, and the neutrino fl ux, respectively 
l|Baron et alJll989t faampp fe Jankall2002|) . 

The fluid velocity v is identified with the local radial 
velocity calculated by the Newtonian code and the metric 
function T is given by 



F 



2m TO v 



(5) 



The velocity-dependent terms were added for a closer 
match with the general relativistic form of t he equation of 
motion I van Riperlll979HBaron et alJll989f) . In the treat- 
ment of neutrino transport general relativistic redshift 
and time dilation effects are included, but for reasons of 
consistency with the Newtonian hydrodynamics part of 
the code the distinction between coordinate and proper 
radius is ignored in the rel ativistic transport equ ations 
(for details, see Sect. 3.7.2 of lRampp fc Jankalfcood) . The 
quality of this approach was ascer tained by a comparison 
with fully relativistic cal culations l|Rampp fc Jankal l2002; 
iLiebendorfer et al-lEio^ . 

In order to calculate the effective relativistic poten- 
tial for multi-dimensional flows we substitute the "spheri- 
cal contribution" 3>(r) to the multi-dimensional Newtonian 
gravitational potential 



$(r,9,ip) 



dr' d9' dip' r' 2 sin 9' 



v 



\r — r 



by the TOV potential <E>tov: 

*cff = $ +$TOV- 



(6) 



(7) 



Here$(r) and$xov are calculated according to Eqs. 
and Pjl. respectively, however with the hydrodynamic 
quantities p, e, P, v and the neutrino quantities E, F, p u 
being replaced by their corresponding angularly averaged 
values. Note that v here refers to the radial component of 
the velocity, only. 
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2.2. Modifications of the TOV potential 

In a recent comparison iLiebendorfer et all l(200fj) found 
that gravity as described by the TOV potential in Eq. (JHJ 
overrates the relativistic effects, because in combination 
with Newtonian kinematics it tends to overestimate the 
infall velocities and to underestimate the flow inertia in 
the pre-shock region. Thus, supposedly via the nonlinear 
dependence of <& c ff on e and P the compactness of the 
proto-neutron star is overestimated, with this tendency 
increasing at later times after core bounce. Consequently, 
the neutrino luminosities and the mean energies of the 
emitted neutrinos are larger than in the corresponding 
relativistic simulation. 

In order to reduce these discrepancies - without sac- 
rificing the simplicity of Newtonian dynamics - we tested 
several modifications of the TOV potential, Eqs. J2J), 
which all act to weaken it. In particular, we studied the 
following variations: 

Case A: In the integrand of Eq. (@J a factor T, Eq. (|SJ) 
is added. Since T < 1 this reduces the gravitational 
TOV mass used in the potential. 

Case B: In Eq. (0J the internal gas energy density and 
the neutrino terms are set to zero, e = E = F = 0, 
which again decreases the gravitational TOV mass. 

Case C: In Eq. Q the internal gas energy is set to zero, 
e = 0, which directly weakens the TOV potential. 

Case D: In the equation for the TOV potential, Eq. J2J), 
TOtov is replaced by ^(rriTOV+m g ). Here a Newtonian 
gravitational mass is defined as m g = m r — with 
the rest mass m r = Air J Q dr' r' 2 p and the mass equiv- 
alent of the binding energy nib = 2n\ dr' r' 2 p $|. As 
m g < totov, the strength of the potential is reduced. 

Case E: Both in the equation for the TOV potential, 
Eq. ©J and the equation for the TOV mass, Eq. QJ, 
we set e = 0. 

Case F: In the equation for the TOV potential, Eq. 10), 
we set r = 1. As T < 1 otherwise, this weakens the 
potential. 

Case G: In the expression for T, Eq. (0, the velocity is 
set to zero, v = 0. Hence, T~ 2 increases in Eq. (J2J. This 
modification is used to also test a potential which is 
even stronger than the unmodified TOV potential. 

In addition to these cases with a modified version of the 
TOV potential, we use the following notations: 

Case N: This denotes the purely Newtonian runs with 
"regular" Newtonian potential. 

Case R: This is the "reference" case with the TOV 
potential as defined by Eq. @. 

Case GR: This case refers to fully relativistic simula- 
tions with either the code CoCoNuT or the Agile- 
BoltzTran neutrino radiation- hydrodynamics code 
of the Oak Ridge-Basel collaboration. 



Note that setting the internal energy density e to zero 
in Case B is unambiguous when a simple EoS is used and 
the particle rest masses are conserved. In general, how- 
ever, particles can be created and destroyed, or bound 
states can be formed (e.g., in pair annihilation processes 
or nuclear reactions, respectively). Then only the sum of 
the rest mass energy and internal energy per nucleon - 
both appear in Eqs. (0 0} only combined in form of the 
"relativistic energy" per unit of mass, (p + e)/p - is well 
defined, but not the individual parts. Therefore there ex- 
ists ambiguity with respect to which contribution to the 
energy is set to zero. In order to assess a possible sensitiv- 
ity of the core collapse results to this ambiguity, we tried 
two different variants of Case B in our Vertex simula- 
tions with microphysical EoS. On the one hand we used 
e = i? = F = 0in Eq. (@J , with e being the internal energy 
density plus an ene r gy no rmalization given by the EoS of 
lLattimer fc Swestvl <jl99 lf> . e = [p + e) — p (m n — A)/m u 
(where to u = 1.66 x 10~ 24 g is the atomic mass unit, m n 
the neutron rest mass, and A = 8.8MeV). On the other 
hand we tested e* = E = F = with e* being e without 
this energy normalization, i.e. e* = (p + e) — pm n /m u . 
Nucleons are then assumed to contribute to the TOV 
mass, Eq. with the vacuum rest mass of the neutron, 
increasing the mass integral and reducing T (Eq. J^J) rela- 
tive to the first case, thus making the effective relativistic 
potential a bit stronger again. In order to compensate for 
this we also set p v = in the TOV potential, Eq. @. 
Both variants are found to yield extremely similar results 
and we therefore will discuss only one of them (the first 
variant) as Case B for the Vertex simulations. 

Ideally, a Newtonian simulation with an effective rela- 
tivistic potential not only yields a solution of the TOV 
structure equations for an equilibrium state (as does 
Case R), but in addition closely reproduces the results 
from a relativistic simulation (Case GR) during a dynamic 
evolution. Applying the modifications of the TOV poten- 
tial listed above, we find that Cases A to D yield improved 
results as compared to Case R, while Cases E to G either 
weaken the potential too much or are very close to Case R. 
These findings are detailed in Sect. ^ and 

2.3. Theoretical motivation 

There are (at least) two basic requirements which ap- 
pear desirable for an effective relativistic potential in a 
Newtonian simulation. Firstly, the far field limit of the 
fully relativistic treatment should be approximated rea- 
sonably well in order to follow the long-term accretion of 
the neutron star and the associated growth of its baryonic 
mass. Secondly, the hydrostatic structure of the neutron 
star should well fit the solution of the TOV equations. 

The second point will be discussed in detail in Sect. 14. ll 
A closer consideration of the first point suggests the mod- 
ified effective relativistic potential of Case A as promis- 
ing, and in fact it turns out to be the most preferable 
choice concerning consistency and quality of the results. 
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The other cases listed in Sect. 12.21 are mostly ad hoc mod- 
ifications of the original effective relativistic potential of 
Eqs. (Case R) with the aim to reduce its strength, 

which was found to overestimate the effects of gravity 
compared to fully relativistic simulations in previous work 
(Liebendorfer et al. 2005). These cases are also discussed 
here for reasons of comparison and completeness. 

In Eq. @ the hydrodynamic quantities (like rest-mass 
density p plus extra terms) are integrated over volume. In 
the Newtonian treatment there is no distinction between 
coordinate volume and local proper volume. Performing 
the integral of Eq. 10} therefore leads to a mass - used 
as the mass which produces the gravitational potential 
in Eqs. and - which is larger than the baryonic 
mass, TOb = 47T J dr' r' 2 p. In particular, it is also larger 
than the gravitational mass in a consistent relativistic 
treatment, which is the volume integral of the total en- 
ergy density and includes the negative gravitational po- 
tential energy of the compact object. The latter reduces 
the gravitational mass relative to the baryonic mass by 
the gravitational binding energy of the star (see, e.g., 
IShapiro fc Teukolskvl [l983, page 125 for a correspond- 
ing discussion). Th erefore, the effect i ve rela tivistic poten- 
tial introduced by iRampp fc Jankal (2002) four Case R, 
Eqs. J2HDl] cannot properly reproduce the far field limit 
of the relativistic case and thus overestimates the effects 
of gravity. This particularly applies to the infall velocities 
of the stellar gas ahead of t he supernova shock, as shown 
in iLiebendorfer et al.N2005h . 

Introducing an extra factor T in the integral of Eq. Q 
for the TOV mass is motivated by the following consider- 
ations (where for reasons of simplicity contributions from 
neutrinos, though important, are neglected and spherical 
symmetry is assumed): In the relativistic treatment the 
total (gravitating) mass of the star is given as 

/>oo 

m g = 4tt / drV 2 (p + e), (8) 
Jo 

where dV = Airdr' r' 2 is the coordinate volume element, 
whereas the baryonic mass is 

/>oo 

m b = 4tt / drV 2 r"V, (9) 
Jo 

with dV' = 47rdrV 2 r~ 1 being the local proper volume 
element. The factor > 1 in the integrand of mb thus 
ensures that mb > m s . The integral for m g can also be 
written as 

/>oo -p />oo 

m g = 47r/ dr' r' 2 - (p + e) = / dVT(p + e). (10) 

Jo r j Q 

Since in Newtonian hydrodynamics no distinction is made 
between coordinate and proper volumes, one may identify 
dV = dV, consistent with the rest of our Newtonian 
code. This leaves the additional factor V in the integrand 
of Eq. (fTO)) . leading to a redefined TOV mass used for 
computing the effective relativistic potential in Case A, 

mTOV ( r )=47r^ dr'r ,2 r (p + e + E + ^ . (11) 



The fact that a factor in the volume integral estab- 
lishes the relation between gravitating mass, Eq. JSJ), and 
baryonic mass, Eq. JjJJl, in the relativistic case suggests 
that the factor V < 1 in Eq. I jlljl might lead to a suitable 
reduction of the overestimated effective potential that re- 
sults when the original TOV mass of Eq. is used in 
Eqs. J5J and (J5J. indeed, a comparison of the integral of 
Eq. (|11|) for large r with the rest-mass energy of a neutron 
star reduced by its binding energy at time t (computed 
from the emitted neutrino energy, f dt' L v {t') with L v 
being the neutrino luminosity) reveals very good agree- 
ment. 

The arguments given above only provide a heuristic 
justification for the manipulation of the TOV potential 
proposed in Case A. A deeper theoretical understanding 
and more rigorous analytical analysis of its consequences 
and implications is certainly desirable, but beyond the 
scope of the present paper. We plan to return to this ques- 
tion in future work. 

3. Numerical methods 

In the following we briefly introduce the two numerical 
codes used for the simulations presented in this work. The 
codes are based on state-of-the-art numerical methods for 
hydrodynamics and neutrino transport coupled to grav- 
ity, and were used to simulate stellar core collapse and 
supernova explosions. 

3.1. Hydrodynamic code without transport 

The purely hydrodynamic simulations are performed with 
the code CoCoNuT developed by IPimmelmeier et"al] 
( 2002a b) with a metric solver based on spe ctral meth- 
ods as described in IPimmelmeier et, alJ l|2005l) . The code 
optionally uses either a Newtonian gravitational potential 
(or an effective relativistic potential, see Sect. |2J), or the 
general relativistic field equations for a curved spacetime 
in the APM 3 + 1-split under the assumption of the con- 
formal flatness condition (CFC) for the three-metric. The 
(Newtonian or general relativistic) hydrodynamic equa- 
tions are consistently formulated in conservation form, 
and are solved by high-resolution shock-capturing schemes 
based upon state-of-the-art Riemann solvers and third- 
order cell-reconstruction procedures. Neutrino transport 
is not included in the code. A simple hybrid ideal gas EoS 
is used that consists of a polytropic contribution describ- 
ing the degenerate electron pressure and (at supranuclear 
densities) the pressure due to repulsive nuclear forces, and 
a thermal contribution which accounts for the heating of 
the matter by shocks: 

P = P p + Pth, (12) 
where 

P p = Kp\ P th = pe th ( 7th - 1), (13) 

and e t h = e — e P - The polytropic specific internal energy e p 
is determined from P p by the ideal gas relation in combi- 
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nation with continuity conditions in the case of a discon- 
tinuous 7. In that case, the polytropic constant K also has 
to be adjusted (for more details, see iDimmelmeier et all 
l2002aH.Ianka et alJll99.1) . 

The CoCoNuT code utilizes Eulerian spherical coor- 
dinates {r, 9, y>}, and thus axially or spherically symmetric 
configurations can be easily simulated. For the core col- 
lapse simulations discussed in Sects. l4~Tl andl5l the finite 
difference grid consists of 200 logarithmically spaced ra- 
dial grid points with a central resolution of 500 m. A small 
part of the grid covers an artificial low-density atmosphere 
extending beyond the core's outer boundary. The spectral 
grid of the metric solver is split into 6 radial domains with 
33 collocation points each. In order to be able to accu- 
rately follow the dynamics, the domain boundaries adap- 
tively contra ct towards the centre dur ing the collapse, as 
described in IDimmelmeier et alJ 1 2005|) . For the migration 
test fSect. ro)) the finite difference grid consists of 250 log- 
arithmically spaced radial grid points, of which 60 and 190 
cover the neutron star and the atmosphere, respectively. 
For the two-dimensional simulations of rotational core col- 
lapse (Sect. |5j) the finite difference and spectral grids are 
extended by 30 and 17 equidistant angular grid points, 
respectively. 

Even when using spectral methods the calculation 
of the spacetime metric is computationally expensive. 
Hence, in the relativistic spherical (rotational) core col- 
lapse simulations the metric is updated only once ev- 
ery 10/1/50 (100/10/50) hydrodynamic time steps be- 
fore/during/after core bounce, and extrapolated in be- 
tween. In the migration test the update is performed every 
10th step during the entire evolution. The numerical ad- 
eq uacy of this procedure is te sted and discussed in detail 
in IDimmelmeier et"al1 ( 2002al) . 



3.2. Hydrodynamic code with neutrino transport 

The Vertex neutrino-hydrodynamics code solves the 
conservation equations of Newtonian hydrodynamics in 
their Eulerian and co nservative form using the PPM-based 
Prometheus code l|Frvxell et all Il990|) . The neutrino 
transport is treated by determining iteratively a solution 
of the coupled set of moments equations and Boltzmann 
equation , achieving closure by a variable Eddington factor 
method l|Rampp fc Jankall2002j) . See the latter reference 
also for details about the coupling between the hydro- 
dynamics and neutrino transport parts of the Vertex 
code. The models presented i n this paper are compute d 
with the equation of state of lLattimer fc Swestvl l)l99lh . 
in agre ement with the choice o f input physics in the 
work of ILiebendorfer et alJ 1 2005|) . where results obtained 
with the Vertex code were compared with those from 
the Newtonian and fully relativistic calculations with the 
Agile-BoltzTran code of the Oak Ridge-Basel collab- 
oration. 

In order to compare the r esults presented he r e with 
those of the calculations of ILiebendorfer et all l|2005h 
we used the same set o f neutrino interact i on rat es as 
picked for Model G15 in ILiebendorfer et alJ l)2005|) . and 
exactly the same parameters for the numerical setup 
(e.g., the grids for hydrodynamics and neutrino trans- 
port). Information abou t this setup can be found in 
ILiebendorfer et all l|2005h . The initial model for our cal- 
culations is the 15Mfn pr ogenitor model "sl5s7b2" from 
IWooslev fc Weaverl ill 9 951) . 

Since solving the neutrino transport problem is compu- 
tationally quite expensive we performed calculations only 
for Cases A, B, and F (as defined in Sect. 12. 2|) with the 
Vertex code. The quality of these results is then com- 
pared to the fully relativistic treatment of the Agile- 
BoltzTran code. 



The quality of the CFC approximation has been tested 



lapse and for neutron star models llCook et alJ 


1996 
2004 


Dimmclmgigr.e^aL 


2002a 


Shibata & Sekieuchi 


Dimmelmeier et al. 


2005; ( 


.^erda-Duran et al.l 


2004). For 



single stellar configurations the CFC approximation is 
very accurate as long as neither the compactness of the 
star is extremely relativistic nor its rotation is too fast and 
differential. Both criteria are well fulfilled during super- 
nova core collapse. Concerning core collapse simulations 
of rapidly and strongly differentially rotating configura- 
tions at very high densities, we recently found that the 
CFC approximation yields excellent agreement with for- 
mulations solvin g the exact spac etime metric even in this 
extreme regime (|Ott et alJ 12(1051) . 

As the three-metric of any spherically symmetric 
spacetime can be written in a conformally flat way, the 
spherical CFC simulations presented in Sect. I4.ll and 14.31 
involve no approximation, i.e., the exact Einstein equa- 
tions are solved in these cases. 



3.3. Hydrodynamics and implementation of the 
effective relativistic potential 

The implementation of an effective relativistic gravita- 
tional potential or of an effective relativistic gravitational 
force as its derivative into existing Newtonian hydrody- 
namics codes is straightforward and does not differ from 
the use of the Newtonian potential or force. 

The equations solved by the two codes used for 
our simulations are described in much detail in pre- 
vious publications (see, e .g.. iRampp fc Jankal 2002; 
Miiller fc Steinmet d Il995t IDimmelmeier et alJ l2002aL 
2005L and references therein). The implementation of the 



source terms of the gravitational potential as discussed in 
Miil ler fc Steinmetzl l|1995l ) is also applied for the handling 
of the effective relativistic potentials investigated in this 
work. The only specific feature here is the mutual depen- 
dence of r and rnxoV) which is either accounted for by a 
rapidly converging iteration or by taking T from the old 
time step in the update of ttitoVj when the changes during 
the time steps are sufficiently small. 



6 



A. Marek et al.: An improved effective potential for supernova simulations 



Since the actual form of the gravitational source term 
is unspecified in the conservation laws of fluid dynamics, 
the Newtonian potential can be replaced by the effective 
relativistic potentials investigated in the current work in 
a technically straightforward way. Solving an equation for 
the internal plus kinetic energy (as in our codes) requires 
a treatment of the gravity source term in this equation 
that is consistent with its implementation in the equation 
of momentum. 

Of course, the effective potential must be investigated 
concerning its consequences for the conservation of mo- 
mentum and energy. Since a potential constructed accord- 
ing to Eqs. J21EJ does not satisfy the Poisson equation, the 
moment um equatio n cannot be cast into a conservation 
form (cf. lShulll992l Part I, Chapter 4). As a consequence, 
the total linear momentum is strictly conserved only when 
certain assumptions about the symmetry of the matter 
distribution are made, for example in the case of spherical 
symmetry or axially symmetric configurations with equa- 
torial symmetry, or when only one octant is modeled in 
the three-dimensional case. In axisymmetric simulations 
the conservation of specific angular momentum is fulfilled 
as well, when using the effective relativistic potential. In 
general, however, a sufficient quality of momentum (and 
angular momentum) conservation has to be verified by in- 
specting the numerical results. 

The long-range nature of gravity prohibits to have 
an equation in pure conservation form for the total en- 
ergy, i.e., for the sum o f internal, kinetic, and gravita- 
tional energy (|Shulll992l Part I, Chapter 4). In contrast 
to the Newtonian case, however, our effective relativis- 
tic potential does also not allow one to derive a conser- 
vation equation for the total energy integrated over all 
space. Monitoring global energy conservation in a simula- 
tion with effective relativistic potential therefore requires 
integration of the gravitational source terms over all cells 
in time. If the local effects of relativistic gravity are con- 
vincingly approximated - as measured by good agreement 
with static solutions of the TOV equation and with fully 
relativistic, dynamical simulations - there is confidence 
that the integrated action of the employed gravitational 
source term approximates well also the global conversion 
between total kinetic, internal, and gravitational energies 
found in a relativistic simulation. 

The recipes for approximating general relativity should 
be applicable equally well in hydrodynamic codes differ- 
ent from our (Eulerian) PPM schemes, provided the ef- 
fects of gravity are consistently treated in the momentum 
and energy equations. The proposed effective potentials 
are intended to yield a good representation of the effects 
of relativistic gravity in particular in the context of stellar 
core collapse and neutron star formation. For our approx- 
imation to work well, the fluid flow should be subrelativis- 
tic. The numerical tests described in the following sections 
show that velocities up to about 20% of the speed of light 
are unproblematic. 



4. Simulations in spherical symmetry 

4.1. Supernova core collapse with simplified equation 
of state 

We first test the quality of the TOV potential and its 
various modifications as an effective relativistic potential 
in the moderately relativistic regime by evolving models 
of supernova core collapse with the simple EoS described 
in Sect. 13.11 Using the CoCoNuT code we are able to 
perform a direct comparison with a Newtonian version 
of the code (employing either the Newtonian, TOV, or 
modified TOV potential) against the relativistic version 
in a curved spacetime. 

The initial models are polytropes, which mimic an iron 
core supported by electron degeneracy pressure, with a 
central density p cl = 10 10 g cm' 3 and EoS parameters 
7i = 4/3 and K = 4.897 x 10 14 (in cgs units). To initi- 
ate the collapse, the initial adiabatic index is reduced to 
7i < 7;. At densities p > p nuc = 2.0 x 10 14 g cm -3 the 
adiabatic index is increased to 72 > 2.5 to approximate 
the stiffening of the EoS. This leads to a rebound of the 
core ( "core bounce" ) . The inner part of the core comes to 
a halt, while a prompt hydrodynamic shock starts to prop- 
agate outward. Since the matter model of the CoCoNuT 
code does not account for neutrino cooling, the compact 
collapse remnant (corresponding to the hot proto-neutron 
star) cannot cool and shrink, but for some models only 
slightly grows in mass due to accretion. These stages of 
the evolution can be seen in the top panel of Fig.^ where 
the time evolution of the central density p c is plotted for 
a typical core collapse model. 

As it is known from numerical models with sophisti- 
cated microphysics that the prompt hydrodynamic shock 
turns into an accretion shock shortly after core bounce, 
we choose the EoS parameters 71 = 1.325, 72 = 3.0, and 
7th = 1-2 in our accretion shock model AS (Fig.^l to re- 
produce such a behaviour. The conversion of the prompt 
shock into an accretion shock at late times is reflected by 
the radial velocity profiles plotted in the bottom panel of 

Fig.m 

When recalculating Model AS with regular Newtonian 
gravity (Case N), the central density remains below the 
one for relativistic gravity (Case GR) during the entire 
evolution (see top panel of Fig. |2J). This behaviour was 
already observed in the rotational core collapse simula- 
tions of iDimmelmeier et alJ l)2002al) . Replacing the regu- 
lar Newtonian potential (Case N) by the TOV potential 
(Case R) yields the opposite effect, i.e., the central den- 
sity is too high throughout the evolution. However, when 
using modifications A or B of the TOV potential, the be- 
haviour of p c agrees very well with that observed in the 
relativistic simulation (Case GR), indicating that the dy- 
namics of the core collapse with this effective relativistic 
potential is similar to that of the relativistic case. 

When applying the same methods to a model with 
a prompt shock which continues to propagate, we ob- 
serve the same qualitative behaviour (see bottom panel 
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Fig. 1. Time evolution of the central density p c (top 
panel) and radial profiles of the velocity v (bottom panel) 
for the accretion shock model AS in relativistic gravity 
(Case GR). The velocity profiles are snapshots long be- 
fore bounce, during bounce shortly before and after max- 
imum density, and long after bounce, respectively. At 
t = 140.0 ms the prompt shock has turned into an ac- 
cretion shock. 
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Fig. 2. Time evolution of the central density p c for the 
accretion shock model AS (top panel) and prompt shock 
model PS (bottom panel). The regular Newtonian simula- 
tion (Case N; dashed line) yields a smaller p c compared to 
the relativistic one (Case GR; solid line). If a TOV poten- 
tial is used (Case R; dashed-dotted line), p c is too high, 
while modifications A (dotted line) and B (long-dashed 
line) of the TOV potential give results close to the rela- 
tivistic one. 



of Fig. [21). For this Model PS the EoS parameters are 
set to 7i = 1.300, 72 = 2.0, and 7th = 1.5, respectively. 
Both during core bounce and at late times (see magnify- 
ing inset) the Newtonian simulations with modifications A 
and B of the TOV potential also yield results which agree 
noticeably better with the relativistic one (Case GR) than 
with the TOV potential (Case R), and obviously also much 
better than with the regular potential (Case N). 

The time plotted in Fig. 121 is t — where t\, is the 
time of bounce for each simulation. We point out that 
the time coordinate t in the relativistic spacetime is the 
time of an observer being at rest at infinity, and thus we 
can directly relate it to the time used in the Newtonian 
simulation (with regular or effective relativistic potential) . 
Moreover, at core bounce the differences between the rel- 
ativistic coordinate time t and the relativistic proper time 
t pc = J a c dt of the fluid element at the center, where a 
is the ADM lapse function, are negligible for the gravita- 
tional fields encountered in our core collapse models. In 
addition, shortly after core bounce the compact collapse 
remnant reaches an approximate hydrostatic equilibrium 
with nearly constant p c . Therefore, for the models con- 
sidered here the time evolution of the central density in 



proper time t p c would look similar to the one in coordinate 
time t. 

The discussion up to now has only been concerned with 
the local quantity p c . While this is a good measure for 
the global collapse dynamics, a direct comparison of, e.g., 
radial profiles at certain evolution times is more powerful. 

In Fig.|3|radial profiles of the density p for the core col- 
lapse model AS are compared at "late" times (t — 150 ms, 
corresponding to about 60 ms after core bounce). As this 
model exhibits no post-bounce mass accretion, shortly af- 
ter core bounce the compact remnant has acquired a new 
equilibrium state with essentially vanishing radial velocity. 
Its density profile can thus be very well described by that 
of a TOV solution constructed with equal central density 
p c and an "effective" one-parameter EoS, P — P(p), ex- 
tracted from the two-parameter EoS, P = P(p, e), defined 
by Eqs. I|12l 1131) . Details of this procedure are given in 
Appendix [X] The top panel of Fig. [3] shows that this be- 
haviour is indeed confirmed by our simulations. The den- 
sity profile obtained from a dynamic simulation with rel- 
ativistic gravity (Case GR) and the one obtained from 
solving the TOV structure equations agree extremely well 
both at supranuclear and subnuclear densities. Only ahead 
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Fig. 3. Comparison of radial profiles of the density p from 
the hydrodynamic evolution of Model AS at late times 
(solid line) and from a solution of the TOV structure equa- 
tions with identical central density (dashed line). In both 
the relativistic simulation (Case GR; top panel) and the 
Newtonian simulation with TOV potential (Case R; bot- 
tom panel) the two curves are almost indistinguishable, 
while the regular Newtonian simulation (Case N; centre 
panel) yields a large mismatch. For comparison, the rela- 
tivistic profile (solid line in the top panel) is also plotted 
in dashed-dotted line style in the lower two panels. 

of the outward propagating shock (off scale in Fig. [3J 
where matter is still falling inward, i.e., where the radial 
velocity is not negligible, a disagreement can be observed. 

When using regular Newtonian gravity (Case N; see 
centre panel of Fig. the situation is different in two 
aspects. First, due to the shallower Newtonian gravita- 
tional potential, the density is significantly smaller in the 
central parts of the compact remnant compared to rela- 
tivistic gravity (Case GR). However, beyond r ~ 10 km 
the density in the Newtonian model exceeds that of the 
relativistic simulation, which is the d ensity crossing e ffect 
observe d and discussed in detail by iDimmelmeier et"afl 
120023). The weaker regular Newtonian potential thus 
yields a remnant which is less compact and dense in the 
centre, but relatively denser in the outer parts. 

Second, when comparing the density profile resulting 
from the hydrodynamic evolution with that obtained by 



r [km] 

Fig. 4. Same as Fig.0 but for radial profiles of the specific 
internal energy e. The step-like shape of the profile arises 
from contributions of the thermal energy eth in various 
strengths. 

solving the TOV structure equations (for the same p c ), a 
clear disagreement can be noticed. As the core reaches a 
final compactness 2M corc /i? coro ~ 0.2 as measured by the 
ratio of core mass to core radius, a hydrodynamic simula- 
tion with regular Newtonian potential cannot be expected 
to yield an equilibrium result which is identical to a solu- 
tion of the (relativistic) TOV structure equations. 

If on the other hand the Newtonian hydrodynamic sim- 
ulation is performed with the TOV potential (Case R; see 
bottom panel of Fig.^J), the final equilibrium state is (as 
in Case GR) a solution of the TOV structure equations. 
However, after core bounce the density profile of Case R 
lies well above that obtained in the relativistic simulation 
(Case GR). Thus Case GR and Case R both exhibit post- 
bounce density stratifications which are consistent with a 
relativistic equilibrium state but otherwise differ notice- 
ably from each other. 

Radial profiles of the specific internal energy (see 
Fig. exhibit a similar behaviour as the density pro- 
files presented in Fig. [3] As expected, both the relativistic 
simulation (Case GR) and the Newtonian simulation with 
TOV potential (Case R) yield final equilibrium profiles 
which closely coincide with the corresponding solution of 
the TOV structure equations (see top and bottom panel, 
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Fig. 5. Comparison of radial profiles of the density p from 
the hydrodynamic evolution of Model AS at late times 
(solid line) and from a solution of the TOV structure equa- 
tions with identical central density (dashed line). In the 
Newtonian simulations with either modification A (top 
panel) or B (bottom panel) of the TOV potential, the two 
curves match well and also agree with the result from the 
relativistic simulation (Case GR; dashed-dotted line). 



respectively). Again, analogous to the density, in Case R 
the internal energy in the compact remnant is significantly 
higher as compared to Case GR. In the Newtonian sim- 
ulation with regular potential (Case N; centre panel) the 
profile deviates clearly from both the corresponding solu- 
tion of the TOV equation with equal central value and the 
profile computed with relativistic gravity (Case GR), and 
does not even show the distinctive steps of that profile. 

The above findings imply that the final equilibrium 
state of a Newtonian core collapse simulation performed 
with the TOV potential (Case R) is consistent with a 
TOV solution of identical central density and specific in- 
ternal energy. However, if relativistic effects are important, 
the density stratification and energy distribution (which 
in turn determines the "effective" EoS as described in 
Appendix EJ will differ from that of a consistently rel- 
ativistic simulation (Case GR). Presumably the combi- 
nation of the relativistic TOV potential (which depends 
nonlinearly on the pressure and the internal energy den- 
sity) with Newtonian kinematics in Case R is responsible 
for this difference in the stratification of density and spe- 
cific internal energy in the compact remnant compared to 
Case GR (as seen in the bottom panels of Figs. |3] and 0J. 
Currently we cannot provide a more detailed analysis for 
understanding the mechanism by which these discrepan- 
cies are caused. We plan to conduct further research on 
this issue. 
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Fig. 6. Same as Fig.|31 but for radial profiles of the specific 
internal energy e. The step-like shape of the profile arises 
from contributions of the thermal energy e t h in various 
strengths. 



As discussed in Sect. 12.21 a modification of a 
Newtonian simulation incorporating the effects of rela- 
tivistic gravity must not only yield equilibrium states 
which (approximately) reproduce a solution of the TOV 
equations, particularly for compact matter configurations. 
More importantly, the density and energy profiles must 
also agree as closely as possible with the profiles obtained 
from the evolution with a relativistic code (which im- 
plies the previous requirement). Both these requirements 
are fulfilled for Newtonian simulations utilizing modifica- 
tion A and B of the TOV potential, as shown for pro- 
files of the density and specific internal energy in the final 
equilibrium state of the core collapse model AS in Figs. [31 
and El respectively. In these cases the profiles agree well 
with the relativistic result (Case GR) not only near the 
centre (as already observed in the time evolution of p c ; see 
top panel of Fig. but also throughout the entire core, 
with the agreement being much better compared to the 
TOV potential (Case R). Accordingly, the solution of the 
TOV structure equations for identical central values yields 
similar profiles as well. 

We conclude that the modifications A and B of the 
TOV potential combined with Newtonian kinematics pro- 
vide a very good approximation of the results obtained 
with relativistic gravity (Case GR) in the strongly dy- 
namic phases during core bounce as well as at post-bounce 
times when the compact remnant acquires equilibrium. 
The same findings hold for the modifications C and D of 
the TOV potential, although we do not present results for 
these cases here. On the other hand, compared to Cases A 
to D, the modifications E to G produce less accurate re- 
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suits. Thus, we consider them inappropriate for an effec- 
tive relativistic potential. 

In passing we note that gauge dependences may play 
an important role when non- invariant quantities (like, e.g., 
density profiles) from simulations with a relativistic code 
are compared against results from either a relativistic code 
using a different formulation of the spacetime metric or 
from a Newtonian code. For instance, the radial coordinate 
r used in the profile plots in Figs. to is the isotropic 
coordinate radius of the ADM-CFC metric used in the 
rela tivistic version of the Co CoNuT code (for details, 
see iDimmelmeier et al.l Eo02a'l. In contrast to this, the 
TOV equations HA.lllA.2jl use the standard Schwarzschild- 
like coordinate radius, which is equivalent to the circum- 
ferential radius. Thus, a coordinate transformation from 
the standard Schwarzschild-like radial coordinate to the 
isotropic radial coordinate (as specified in Appendix [5| 
is applied to all profiles of quantities obtained with the 
TOV solution or Newtonian simulations in those figures. 
Otherwise, the difference in the two radial coordinates of 
up to ~ 20% would lead to a noticeable distortion of the 
radial profiles. We also point out that we neglect any in- 
fluence of the difference between Newtonian and relativis- 
tic coordinate time (or proper time) when comparing the 
profiles at the same coordinate times, as the compact rem- 
nants are in equilibrium at "late" times (about 60 ms after 
core bounce). 



4.2. Supernova core collapse with microphysical 
equation of state and neutrino transport 

Next we explore the moderately relativistic regime of 
stellar core collap s e wit h the microphysical EoS of 
lLattimer fc Swestvl l|l99l|) and neutrino transport. We 
simulated the collapse and the post-bounce evolution of 
the progenitor model sl5s7b2 with the Vertex code as 
detailed in Sect. 18.21 The calculations were performed 
using the TOV potential given in Eq. © (Model V-R, 
which i s identical with the Vert ex calculation of Model 
G15 in ILiebendorfer et al.ll2005^ . and we also tested the 
modifications A, B, and F of the TOV potential (Models 
V-A, V-B, and V-F, respectively). For comparison, we 
ref er in the following discus sion also to the calculation 
of ILiebendorfer et all l|2005h with the fully relativistic 
Agile-BoltzTran code (Case GR, Model AB-GR). 

Fig. shows the central density as a function of time 
for the collapse (left panel) and for the subsequent post- 
bounce phase (right panel). The "central" density is the 
density value at the centre of the innermost grid zone of 
the AB-GR simulation. Because of a different numerical 
resolution it was necessary to interpolate the Vertex re- 
sults to this radial position. During the collapse only mi- 
nor differences between the relativistic calculation (bold 
solid line) and the calculations with the Vertex code are 
visible. Note that the trajectories from the Vertex code, 
with the modifications A, B, and F of the TOV potential 




-100 
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Fig. 7. Time evolution of the central density p c for Model 
AB-GR (bold solid line), V-A (solid line), V-B (dashed- 
dotted line), V-F (dashed line), and V-R (dotted line). 
The left panel shows the collapse phase (note that here all 
models with the Vertex code lie on top of each other) , 
while the right panel shows the post-bounce evolution. 
Note the different axis scales in both panels. 



as well as with the TOV potential (case R), lie on top of 
each other. 

We can thus infer that the differences between the 
modifications A, B, and F of the TOV potential are unim- 
portant during the collapse phase. Furthermore, we can 
conclude from the good agreement of the general relativis- 
tic calculation and the Vertex calculations that the TOV 
potential works well during the collapse phase. However, 
after core bounce this potential overestimates the com- 
pactness of the forming neutron star (just as in the purely 
hydrodynamic simulations in Sect. I4.l|l . and therefore the 
density trajectories of Model AB-GR and Model V-R di- 
verge (see Fig.[7J|. At 250 ms after the shock formation the 
central density in Model V-R is about 20% higher than the 
one in the relativistic calculation. At this time the mod- 
ifications A and B of the TOV potential give a central 
density only about 2% higher than Model AB-GR, and 
the absolute difference stays practically constant during 
the entire post-bounce evolution. This implies that both 
modifications yield very good quantitative agreement with 
the general relativistic treatment. In contrast, in Model 
V-F the central density after bounce is lower than the rel- 
ativistic result of Model AB-GR. This indicates a strong 
underestimation of the depth of the gravitational potential 
in Case F, where T = 1 in the integrand of Eq. 

Since the central densities suggest that differences 
between a fully relativistic calculation and Newtonian 
simulations with effective relativistic potential be- 
come significant only a fter shock formation (see also 
ILiebendorfer et aLll2005|) . we discuss the implications of 
our potential modifications in the following only during 
the post-bounce evolution. 
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Fig. 8. Time evolution of the shock position r s for simu- 
lations with the Vertex code using various modifications 
of the TOV potential compared to the general relativistic 
result from the Agile-BoltzTran code (Model AB-GR, 
thick solid line). 

Fig. |Hl shows the shock positions as functions of time. 
Both Case A (thin solid line) and Case B (dashed-dotted 
line) reveal the desirable trend of a closer match with the 
general relativistic calculation (thick solid line) than seen 
for Model V-R, which gives a shock radius that is too 
small, and Model V-F, where the shock is too far out at 
all times. In particular, Model V-A reveals excellent agree- 
ment with Model AB-GR. The only major difference is 
visible between 170 ms and about 230 ms after shock for- 
mation when the shock transiently expands in the Vertex 
calculation. This behaviour is generic for the Vertex re- 
sults and independent of the choice of the gravitational po- 
tential. In the Agile-BoltzTran run the transient shock 
expansion is much less pronounced and also a bit delayed 
relative to the Vertex feature (it is visible as a deceler- 
ation of the shock retraction between about 200 ms and 
250 ms). This difference, however, is not caused by gen- 
eral relativistic effects but is a consequence of a different 
numerical tracking of the time evolution of interfaces be- 
tween composition layers in the collapsing stellar core (for 
more details about the numerics and a discu ssion of the 
involved physics, see lLiebendorfer et al.ll2005)) . It is there- 
fore irrelevant for our present comparison of approxima- 
tions to general relativity. A good choice for the effective 
relativistic potential (like Case A) should just ensure that 
the corresponding shock trajectory converges again with 
the relativistic result (Case GR) after the transient period 
of shock expansion. 

The time evolution of the central density or the shock 
position, however, is not the only relevant criterion for as- 
sessing the quality of approximations to general relativity, 
as already discussed in Sect. 14.11 A good approximation 
does not only require good agreement for particular time- 
dependent quantities (like, e.g., the central density), but 
also requires that the radial structure of the models re- 



produces the relativistic case as well as possible at any 
time. 

In the left panels of Fig. [§| we show such profiles of 
the density and velocity (top panel) and of the entropy 
and electron fraction Y e , i.e., the electron-to-baryon ratio 
(bottom panel), for Models AB-GR and V-A at a time 
of 250 ms after bounce, when the discrepancy between 
relativist ic and approximative trea tment was found to be 
largest in lLiebendorfer et alJ l)2005|) . Fig.^can be directly 
compared with Fig. 12 in the latter reference. Obviously 
Model V-A fits the density profile of the general relativis- 
tic calculation (Model AB-GR) extremely well at all radii. 
Furthermore, both the velocity ahead of the shock front 
and behind it are in extremely good agreement between 
the two models (the differences at the shock jump have 
probably a numerical reason associated with the differ- 
ent handling of shock discontinuities in both codes). It is 
an important result that with this modification A of the 
TOV potential one is able to approximate the kinematics 
of the relativistic run with astonishingly good quality in a 
Newtonian calculation (at least in supernova simulations 
when the velocities do not become highly relativistic). In 
contrast, in Vertex runs with the original TOV potential 
(Case R) the pre- shock velocities were foun d to be signif- 
icantly too large ^Liebendorfer et al"]l2005[) . which is due 
to the overestimated strength of gravity in the far field 
limit as compared to the relativistic calculations (see the 
discussion in Sect. 12.31) . 

Also the entropy and Y e profiles (bottom left panel 
of Fig. reveal a similarly excellent agreement between 
Models V-A and AB-GR. The minor entropy differences 
ahead of the shock are associated with a slightly differ- 
ent description of the microphysics (nuclear burning and 
eq uation of state) in the in fall region (for details we refer 
to iLiebendorfer et alJl2005T) . 

Not only the radial structure of the forming neutron 
star in all relevant quantities is well reproduced, but also 
the neutrino transport results of the relativistic calcula- 
tion and of the approximative description of Case A are 
in nearly perfect agreement. Corresponding radial profiles 
of the luminosities an d root mean square energ ies - both 
as defined in Sect. 4 of lLiebendorfer et alJ l|2005f) - for elec- 
tron neutrinos, v ei electron antineutrinos, v e , and heavy- 
lepton neutrinos 1 , v x , are displayed in the right panels of 
Fig. H The results for Models V-A and AB-GR for all 
neutrino flavours share their characteristic features, and 
in particular agree in the radial positions where the differ- 
ent luminosities start to rise. While the v x luminosities are 
nearly indistinguishable below the shock, the jump at the 
shock is slightly higher for the Vertex run and reflects 
the larger effects due to observer motion, e.g., Doppler 
blueshift and angular aberration, for an observer comoving 
with the rapidly infalling stellar fluid ahead of the shock. 
The offset between results of Models V-A and AB-GR de- 



1 Since the transport of muon and tau neutrinos and antineu- 
trinos differs only in minor details we treat all heavy-lepton 
neutrinos identically in the Vertex simulations. 
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Fig. 9. Left: Radial profiles of the velocity v (dashed lines) and density p (solid lines) for Model V-A (thin) and the 
relativistic Model AB-GR (bold) at a time of 250 ms after shock formation (top panel), as well as radial profiles of the 
entropy s (da shed lines) and of the elec tron fraction Y e (solid lines) for the same models and the same time (bottom 
panel). As in ILiebendorfer et al.l l)2005|) . r is the circumferential radius in case of the relativistic results. Note the 
different vertical axes on both sides of the two panels. Right: Radial profiles of the luminosities L of electron neutrinos 
(solid lines), electron anti- neutrinos (dotted lines), and heavy- lepton neutrinos (dashed lines) for Models V-A (thin) 
and AB-GR (bold) at a time of 250 ms after shock formation (top panel), as well as radial profiles of the root mean 
square energies (c 1 /)rms for the number densities of v ei D ei and heavy-lepton neutrinos for Models V-A and AB-GR 
(bottom panel). The labeling is the same as in the panel above, and all neutrino quantities are given for a comoving 
observer. 



creases at larger radii where the infall yelocities are lower. 
This discrepancy was not discovered by Liebendorfer et al. 
|)2005[K because there the agreement of the radial struc- 
ture for both investigated models was generally found to 
be poorer than in the present work. 

General relativistic effects are unlikely as an explana- 
tion, beca use they are very small a round the shock (see 
Fig. 13 in ILiebendorfer et al.ll2005j) . A detailed analysis 
reveals that both codes produce internally consistent re- 
sults, conserving to good precision the luminosity through 
the shock for an observer at rest and showing the ex- 
pected and physically correct behavior in the limit of 



large radii. Most of the observed difference (which has 
no mentionable significance for supernova modeling) could 
be traced back to the fact that Vertex achieves only or- 
der (v/c) accuracy, whereas Agile-BoltzTran produces 
the full relativistic result including higher orders in (v/c). 
Corresponding effects become noticeable when v/c <; 0.1. 
The mean neutrino energies are hardly affected by this 
difference (Fig. [^1 bottom right panel). In case of the v e 
and v e luminosities the Vertex run yields roughly 10% 
lower values outside of the corresponding neutrino spheres 
(i.e., between about 50 km and 90 km), but values much 
closer to those from the Agile-BoltzTran calculation 
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Fig. 10. Luminosities as functions of post-bounce time 
for different cases computed with the Vertex code and 
for Model AB-GR. The top panel shows the results for 
heavy-lepton neutrinos, the centre panel those for the elec- 
tron antineutrinos, and the bottom panel the results for 
electron neutrinos. The panels on the left magnify the 
early post-bounce phase. All luminosities are given for an 
observer comoving with the stellar fluid at a radius of 
500 km. Note the different scales of the vertical axes. 



ahead of the shock. Since the neutrinospheric emission 
of v e and v e is strongly affected by the mass accretion 
rate of the nascent neutron star and the corresponding 
accretion luminosity (which both seem to have the ten- 
dency of being slightly higher in Model AB-GR), we re- 
frain from ascribing the different magnitude of the v e and 
v e luminosities only to the treatment of relativistic effects. 
Although such a connection cannot be excluded, the lu- 
minosity differences might (partly) also be a consequence 
of the different accretion histories in Models AB-GR and 
V-A, which manifest themselves in the shock trajectories 
(Fig. |HJ) and are attributable to the different handling of 
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Fig. 11. Radial profiles of the velocity v (dashed lines) 
and density p (solid lines) for Models V-A (bold lines), 
V-B (thin), and V-R (medium) at a time of 250 ms after 
shock formation (top panel) . Radial profiles of the entropy 
s (dashed lines) and of the electron neutrino luminosity 
L Ue (solid lines) for the same models and the same time 
(bottom panel). The luminosity is given for an observer 
comoving with the stellar fluid. 

the microphysi cs and computational grid in the infall layer 
(see above and lLiebendorfer et al.ll2005^ . This interpreta- 
tion seems to be supported by the time evolution of the 
neutrino luminosities plotted in Fig. 1101 The accretion 
bump in the v e and v e luminosities which follows after 
the prompt v e burst is stretched in time in case of Model 
AB-GR, indicating the delay of mass infall at higher rates 
relative to all Vertex simulations. Note that the neutrino 
emission reacts with a time lag of some 10 ms (correspond- 
ing to the cooling timescale of the accretion layer on the 
neutron star) to variations of the mass accretion rate. 

Moreover, Fig. 1101 shows that our variations of the ef- 
fective relativistic potential in the Vertex models have 
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little influence on the prompt burst of v e at shock break- 
out. But subsequently the overestimated compactness of 
the proto-neutron star in Model V-R, which causes the 
faster contraction of the stalled shock after maximal ex- 
pansion (Fig. |SJ) , also leads to higher neutrino luminosi- 
ties during the accretion phase. Consistent with the shock 
trajectories, Model V-A yields the closest match with the 
general relativistic run of Model AB-GR also for the neu- 
trino luminosities. It is very satisfactory that the results 
(shock radius r s as well as the neutrino luminosities L) 
from both simulations reveal convergence at later times 
when the period of massive post-bounce accretion comes 
to an end. 

In Fig. ^2 w e present the radial structure at 250 ms 
after bounce for the Vertex simulations with the modi- 
fications A and B of the TOV potential, compared to the 
results with the T OV potential (Case R) w hich was al- 
ready discussed in ILiebendorfer et al.1 l|2005|) . Note that 
because of the excellent agreement seen in Fig. [51 Model 
V-A (Case A) can also be considered as a representation 
of the fully relativistic run of Model AB-GR. Models V-A 
and V-B show results of similar quality. The little offset 
of the shock position (which is causally linked to the dif- 
ferences in all profiles) might suggest that Model V-B is 
slightly inferior to Model V-A in approximating relativ- 
ity. This conclusion could also be drawn from the post- 
bounce luminosities in Fig. 1101 However, caution seems to 
be advisable with such an interpretation, being aware of 
the uncertainties in the accretion phase and infall layer 
discussed above, and in view of the fact that the central 
densities (Fig. 2} and radial density profiles fFig. 111(1 agree 
well. Moreover, the quality of the agreement at "very late" 
times cannot be judged, because no information is avail- 
able for the behaviour of Model AB-GR after 250 ms post 
bounce, a time when the settling of the shock radius and 
luminosities to their post-accretion levels seems not yet 
over in this model (Figs. El EJ. The TOV potential of 
Case R clearly produces too large infall velocities ahead 
of the shock (and therefore does not agree well with the 
kinematics of the relativistic calculation), overestimates 
the compactness of the forming neutron star, and thus 
produces too high neutrino luminosities during the sim- 
ulated period of evoluti on (for a detailed discussion, see 
ILiebendorfer et a"fll2005f) . Cases A and B clearly perform 
better and must be considered as significant improvements 
for use in Newtonian simulations with an effective rel- 
ativistic potential as approximations to fully relativistic 
calculations. 

4.3. Neutron stars in equilibrium - Migration test 

The supernova core collapse simulations presented in 
Sect. 14.11 and 14.21 are limited to a maximum core com- 
pactness 2M coro /i? corc ~ 0.2. In order to extend the as- 
sessment of the quality of a Newtonian simulation with 
an effective relativistic potential to a more strongly rela- 
tivistic regime, we construct spherical equilibrium models 



of neutron stars using a simple poly tropic EoS, P = Kp 7 , 
with 7 = 2 and K = 1.455 x 10 5 (in cgs units). The TOV 
equations are solved on a very fine grid with 10 5 equidis- 
tantly spaced radial grid points using a second-order ac- 
curate Runge-Kutta integration scheme. By varying the 
value of the central density from p c = 1.0 x 10 14 g cm -3 
to 100.0 x 10 14 g cm -3 , we cover the entire range from 
weak to strong relativistic gravity. The density varia- 
tion corresponds to compactness parameters ranging from 
2M g /R ~ 0.02 to 0.5, where M g is the total gravitational 
mass of the neutron star and R its radius. We again point 
out that in spherical symmetry the assumption of CFC 
constitutes no approximation. 

The mass-density curve of the corresponding TOV 
neutron star models (Fig. El top panel) exhibits the well- 
known form consisting of a stable branch at low densi- 
ties and an unstable branch at high densities, where the 
gravitational mass M g increases or decreases with p c , re- 
spectively. The two branches meet at the maximum mass. 
Configurations on the unstable branch are unstable to 
small perturbations and either collapse to a black hole 
or migrate to a configuration on the stable branch with 
the same gravitational mass but a lower central density 
Pet < Pci (as indicated by the arrow for an exemplary 
configuration in the top panel of Fig. IT2"|) . 

While the gravitational mass M g and the rest mass 
M r are limited in relativistic gravity (Case GR), M g ~ 
M r ~ 1.6M Q , regular Newtonian gravity (Case N) allows 
for solutions with arbitrarily high rest mass M r , which are 
all stable against (small) perturbations. If an (unphysi- 
cal) equivalent to the relativistic gravitational mass is in- 
troduced in Newtonian gravity, M g = M r — Mb, where 
Mb = |i?b| is the mass associated with the binding en- 
ergy of the self-gravitating object, we find negative values 
at large densities, where the binding energy becomes so 
large that M b > M r . 

On the other hand, when using the modifications A 
and B of the TOV potential in the TOV equation for the 
neutron star model sequence, the corresponding curves for 
the TOV gravitational mass M g in Fig.ll2lagree very well 
with the relativistic curve (Case GR) even in the high den- 
sity regime, and they show the same qualitative behaviour 
with an increasing and a decreasing branch. The respec- 
tive mass maxima are located close to the relativistic one 
for both modifications of the TOV potential. 

Note that for constructing equilibrium configurations 
using the TOV potential (Case R) is identical to solv- 
ing the fully relativistic TOV equations (Case GR), and 
thus the results are identical. Therefore, this test is use- 
less for predicting the quality of the TOV potential in the 
context of a dynamic evolution. This is different when a 
neutron star model from the unstable branch is allowed to 
evolve in time using a dynamic code. Then the (small but 
nonzero) finite difference discretisation errors of the nu- 
merical scheme act as perturbations of the unstable equi- 
librium and excite small oscillations. Depending on the nu- 
merical algorithms used in the code and the resolution of 
the grid, the neutron star either contracts to higher densi- 
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Fig. 12. Mass-density diagram for spherical equilibrium 
neutron star models (top panel) and time evolution of the 
central density for the migration test from the unstable 
to the stable branch (bottom panel). In both graphs the 
behaviour of the Newtonian models computed with the 
modification A (dotted line) and B (long-dashed line) of 
the TOV potential is close to that in the relativistic case 
(Case GR; solid line). The Newtonian equivalent of the 
gravitational mass M g (dashed line) as well as the rest 
mass M r (dashed-dotted line) calculated using a regular 
Newtonian potential (Case N) deviate strongly from the 
relativistic solution (top panel). While p c oscillates around 
its initial value p cl for Case N (dashed line), in Case R 
(dashed-dotted line) the neutron star disperses to a state 
with p c <C p c f (bottom panel). 



ties (and ultimately to a black hole for a simple polytropic 
EoS) or approaches a new equilibrium state on the sta- 
ble branch after a strongly dynamical, nonlinear evolution 
phase. This neutron star migration scenario is a standard 
test of the relativistic regime with strongly dynamical evo- 
lution for numeri cal hydrodynamic codes with relativi stic 
gravity (see, e.g.. lFont et "aLll2002t iBaiotti et alJ^OO^ - 



The time evolution of the central density is depicted in 
the bottom panel of Fig.^]for such a migration to the sta- 
ble branch, where the central density p c drops from p c l = 
49.3 x 10 14 g cm~ 3 to p ci = 8.8 x 10 14 g cm- 3 . Shown 
is the evolution for a relativistic simulation (Case GR) 
as well as for Newtonian simulations utilizing a regu- 
lar Newtonian potential (Case N), the TOV potential 
(Case R), and its modifications A and B. For these sim- 
ulations we again use the CoCoNuT code, whose imple- 
mentation of the hybrid EoS, Eqs. (|12l I13jl allows one to 
suppress the contribution of the thermal pressure Pth, in 
which case the hybrid EoS reduces to a polytropic one. 

As expected and as observed in simulations with simi- 
lar relativistic codes, in relativistic gravity (Case GR) the 
neutron star initially expands rapidly and its central den- 
sity decreases strongly. It undershoots the stable equilib- 
rium at p c f and experiences several ring-down oscillations 
of decreasing amplitude until p c approaches p c t at late 
times. The modifications A and B of the TOV potential 
both have initial (final) states on the unstable (stable) 
branch of the mass-density curve which are very close to 
the relativistic one (Case GR) . Therefore the evolution of 
p c (starting again with the same central density p cl ) is 
similar to the relativistic shown in Fig. ^] The 

disagreement in oscillation amplitudes and periods can be 
attributed to the differences in the Newtonian and rela- 
tivistic kinematics. When comparing the results of migra- 
tion tests using variou s relativistic hydrodynamic codes 
l|Stergioulas et al1l2003fl . we observe that even small vari- 
ations in the numerical approach, e.g., a different coordi- 
nate choice or grid setup, or differences in the treatment 
of the artificial atmosphere, can have a strong impact on 
the detailed evolution of p c from the initial to the final 
equilibrium state. Consequently, we consider the results 
obtained with the modifications A or B of the TOV po- 
tential in an otherwise Newtonian formulation to be a very 
good approximation of the relativistic one (Case GR). 

Using regular Newtonian gravity (Case N), however, 
yields a totally different behaviour, because there exist 
no unstable configurations in that case. Thus, discreti- 
sation errors cannot trigger a migration, and the star 
simply oscillates around its initial equilibrium state with 
p c (t) ~ pci- On the other hand, using the TOV potential 
(Case R) leads to a rapid dispersion (without ring-down 
oscillations) of the neutron star with a final central density 
Pc *C p c f, which is in strong qualitative disagreement with 
the behaviour in relativistic gravity (Case GR). We thus 
conclude that in the demanding test case of a migrating 
spherical neutron star model, both the regular Newtonian 
potential (Case N) and the TOV potential (Case R) fail, 
while the modifications A and B of the TOV potential re- 
produce the relativistic result very well, both qualitatively 
and also quantitatively. 

5. Multi-dimensional simulations with rotation 

For the multi-dimensional calculations we only use the 
CoCoNuT code, since parameter studies with MudBath 
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(the two-dimensional version of Vertex) are compu- 
tationally too expensive when neutrino transport is in- 
cluded. Moreover, currently there exists no other fully rel- 
ativistic two-dimensional code including neutrino trans- 
port to compare our results with. The following discussion 
hence focuses on purely hydrodynamic simulations using 
the simple matter model described in Sect. 13.11 

Motivated by the promising results in spherical sym- 
metry, we now apply the modified effective relativistic po- 
tentials according to Eq. J7J to simulations of rotational 
core collapse in axisymmetry. The initial models of the 
pre-collapse core are described by the same EoS parame- 
ters as the ones in Sect. 14. ll with p c \ = 10 10 g cm -3 , but 
now rotate at various rotation rates with different degrees 
of differential rotation. To construct these rotating equilib- 
rium config urations, we us e Hachisu's self-consistent field 
method (see lKomatsu et al.ll989|) . The initial pressure re- 
duction to trigger the collapse and the stiffening of the EoS 
follows the prescription given in Sect. POl 

For these multi-dimensional tests we choose three rota- 
tional core collapse models as representative cases. Their 
evolution is influenced by the effects of rotation at vari- 
ous degree s (for details about the sp ecifications of these 
models, see lDimmelmeier et all2f)02al) . The core of Model 
A1B3G1 is an almost uniform rotator with a moderate 
rotation rate and collapses nearly spherically, forming a 
compact remnant immediately after core bounce. Model 
A4B1G2 is a very differentially rotating model, which de- 
velops considerable deviations from spherical symmetry 
during its evolution, but otherwise shares the qualitative 
bounce behavior of the previous model A1B3G1. Finally, 
in the multiple bounce model A2B4G1 the increase of the 
centrifugal forces during the contraction phase causes a 
series of subsequent multiple bounces at subnuclear den- 
sities. The time evolution of the central density for these 
three models computed relativistically (Case GR) is shown 
with solid lines in the top, centre, and bottom panels of 
Fig. El respectively. 

In a Newtonian simulation with regular Newtonian 
gravity (Case N) all three models exhibit several consecu- 
tive bounces (at subnuclear densities for Models A1B3G1 
and A2B4G1, and just above nuclear matter densities for 
Model A4B1G2), with the central density being consid- 
erably lower during the entire evolution (Fig. I13JI . Such 
a qualitative change of the collapse dynamics between 
a Newtonian and a general relativistic simulation of the 
same model was also observed in the comprehensive pa- 
rameter study of rotationa l core collapse performed by 
iDimmelmeier et alJ l)2f)f)2bh . 

When applying the TOV potential (Case R) or its 
modifications A and B to Model AIB3GI (top panel of 
Fig. 113(1 the outcome is similar to the situation in spher- 
ical symmetry shown in Fig. [3 Case R yields a central 
density noticeably higher than the one in the relativistic 
simulation (Case GR). For Case A the result agrees very 
well with the relativistic one, while for Case B the central 
density is too small. 
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Fig. 13. Time evolution of the central density p c for 
the rotating core collapse models A1B3G1 (top panel), 
A4B1G2 (centre panel), and A2B4G1 (bottom panel), re- 
spectively. With increasing influence of rotation on the col- 
lapse dynamics, the results of the Newtonian simulations 
with the modifications A (dotted line) and B (long-dashed 
line) of the TOV potential and eventually also with the 
TOV potential (Case R; dashed-dotted line) move away 
from the relativistic result (Case GR; solid line). The 
central density obtained with regular Newtonian gravity 
(Case N; dashed line) is lowest for all three models. 

As the influence of rotation increases, the ordering of 
the curves for Cases R, A, and B (from higher to lower p c ) 
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remains unaltered. However, all curves are shifted to lower 
densities (Fig.H3J. For Model A4B1G2 the result from the 
Newtonian simulation with the TOV potential (Case R) 
matches the relativistic result (Case GR) best, which is a 
coincidence associated with the particular rotational state 
of this model. For Model A2B4G1, where the influence of 
rotation is stronger, the results obtained for Cases A, B, 
and R differ considerably from the relativistic results. 

A behaviour consistent with that of these three rep- 
resentative models is observed for other rotational core 
collapse models as well. Thus, in rotational core col- 
lapse, even a Newtonian simulation with the TOV po- 
tential (Case R) yields lower central densities than the 
corresponding relativistic simulation (Case GR) beyond a 
certain rotation rate (which depends both on the initial 
model and on the EoS parameters). 

We therefore conclude that for rotational core collapse 
the use of neither the modifications A and B of the TOV 
potential nor the (stronger) TOV potential itself (Case R) 
can closely reproduce the collapse behaviour of a relativis- 
tic code for all degrees of rotation. Hence, in a Newtonian 
code either the effective relativistic potential has to in- 
clude the effect of rotation, or the collapse kinematics (and 
thus the hydrodynamic equations) have to be modified. 
However, in the latter case the advantage of a Newtonian 
approach over a fully relativistic formulation in terms of 
simplicity and computational speed may be lost. 

On the other hand, we find that the qualitative dynam- 
ical behaviour of relativistic rotational core collapse mod- 
els (Case GR), particularly their collapse type (regular 
collapse, m ultiple bounce collapse, o r rapid collapse as in- 
troduced bv lZwerger fc MiillerllOT^ is correctly captured 
in a Newtonian simulation for a wide range of parameter 
space, if the TOV potential (Case R) or its modified ver- 
sions (Cases A to D) are used. Moreover, the correspond- 
ing results agree much better with the relativistic result 
than that obtained from a Newtonian simulation with reg- 
ular Newtonian potential (Case N). Hence, the usage of 
the effective relativistic potentials is an improvement of 
Newtonian simulations also in the case of rotational core 
collapse. 

6. Summary and Conclusions 

We investigated different modifications of the TOV poten- 
tial used as effective relativistic potentials in Newtonian 
hydrodynamics with the aim to test their quality of ap- 
proximating relativistic effects in stellar core collapse and 
supernova simulations. This work was motivated by a re- 
cent comparison of the neutrino radiation-hydrodynamics 
codes Agile-BoltzTran and Vertex used by the Oak 
Ridge-Base l collaboration and b y the Garching group, re- 
spectively IjLiebendorfer et alJ I2005I) . While the former 
code is a fully relativistic implementation of neutrino 
transport and hydrodynamics in spherical symmetry, the 
latter code employs an approximative relativistic descrip- 
tion by treating the self-gravity of the stellar fluid with 
an effective relativistic potential tailored from a compar- 



ison of the Newtonian and relativistic equations of mo- 
tion. Otherwise Vertex solves the Newtonian equations 
of hydrodynamics, ignoring the ef fects of relativistic kine - 
matics and space-ti me curvature llRampp fc ,Tankall2002h . 
The comparison bv lLiebendorfer et alJ l|200l4l shows very 
good agreement of the results from both codes for the 
phases of stellar core collapse, bounce and prompt shock 
propagation. However, it also revealed a moderate but 
gradually increasing overestimation of the compactness of 
the nascent neutron star during the subsequent accretion 
phase, accompanied by an overestimation of the luminosi- 
ties and mean energies of the emitted neutrinos. 

These systematic differences suggested that an im- 
provement might be possible by changing the effective 
relativistic potential such that the strength of the grav- 
itational field is reduced relative to the TOV potential. 
To this end we considered a variety of modifications 
(Cases A-F; see Sect. of the TOV potential (which 
defines our reference Case R), and compared correspond- 
ing hydrodynamic as well as hydrostatic solutions with 
regular Newtonian (Case N) and fully relativistic calcula- 
tions (Case GR). While some modifications did not pro- 
vide any significant advantage over the TOV potential, 
several options for the potential were identified which al- 
low for a much better reproduction of relativistic results 
with Newtonian hydrodynamics. In particular, reducing 
the TOV mass, Eq. that appears in the expression 
of the TOV potential, turned out to be very successful in 
this respect. 

For example, in Case A we introduced an additional 
metric factor T < 1, Eq. in the integrand of the TOV 
mass, and in Case B we ignored the internal energy (and 
all terms depending on neutrinos) in the TOV mass. In 
both cases not only characteristic model parameters like 
the central density and shock radius were in excellent 
agreement with relativistic simulations during all stages 
of the evolution including the long-term accretion phase 
of the collapsing core after bounce. Also the radial profiles 
of the hydrostatic central core as obtained in the hydrody- 
namic simulations were found to be very good approxima- 
tions of the profiles from the fully relativistic treatment. 
Consequently, they are nearly consistent with solutions 
of the TOV equations for equal central density. Migration 
tests of extremely relativistic, ultradense, hydrostatic con- 
figurations from the unstable to the stable branch of the 
mass-density relation confirmed the good qualitative and 
quantitative performance of the approximation even for 
such a demanding test problem with a strongly dynamic 
transition from the initial to the final state. This suggests 
that the use of a modified TOV potential in an other- 
wise Newtonian hydrodynamics code also ensures a good 
reproduction of relativistic kinematics. 

In addition to these purely hydrodynamic simulations 
with the CoCoNuT code, we performed runs including 
neutrino transport and compared them to relativistic cal- 
culations with the Agile-BoltzTran code of the Oak 
Ridge- Basel collaboration (cf. t he simulations for Model 
G15 in lLiebendorfer et alJl200iA . Also in this application, 
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the modifications A and B of the TOV potential pro- 
duced a major improvement relative to the original TOV 
potential (Case R), for which results with the Vertex 
code were already presented in lLiebendorfer et alJ l|2005jh 
Good agreement between the relativistic and approxima- 
tive treatments was seen in the transport quantities as 
well. The Vertex code performs well despite the fact 
that the approximative description takes into account only 
relativistic redshift in the transport but - for reasons of 
consistency with the Newtonian hydrodynamics code - 
ignores the difference between proper radius and coordi- 
nate radius. Remaining (minor) quantitative differences 
(below about 10% for neutrino data and even smaller for 
hydrodynamic quantities) are probably only partly asso- 
ciated with relativistic effects. We suspect that a signifi- 
cant contribution to these differences originates from small 
(but unavoidable because linked to the specifics of the 
employed numerical methods) discrepancies in the infall 
layer ahead of the shock as described by the Vertex and 
Agile- BoltzTran codes. 

The results discussed in this paper therefore demon- 
strate that a very good, simple, and computationally ef- 
ficient approximation to a relativistic treatment of stel- 
lar core collapse and neutron star formation in spherical 
symmetry can be achieved by modeling the self-gravity of 
the stellar plasma with an effective relativistic potential 
based on the TOV potential as suggested in this work. 
Two-dimensional axisymmetric collapse calculations (us- 
ing a simple equation of state and ignoring neutrino trans- 
port) for rotating stellar cores with a wide range of condi- 
tions showed that even in this case the TOV potential and 
its modifications reproduce the characteristics of the rel- 
ativistic collapse dynamics quantitatively well for not too 
rapid rotation. There is still good qualitative agreement 
when the rotation becomes fast, in contrast to Newtonian 
simulations with regular Newtonian potential (Case N), 
which mostly fail even qualitatively. Multi-dimensional 
simulations with an effective relativistic potential (which 
in general does not satisfy the Poisson equation), however, 
fulfill strict momentum conservation only for special cases 
of symmetry (for details, see Sect. 13.3)) . 

Although our Cases A and B yield results of similar 
quality for all the test problems considered here and nei- 
ther is preferred when a simple equation of state is used, 
the situation is more in favour of Case A when micro- 
physics is taken into account in the equation of state. 
When particles can annihilate or nuclear interactions take 
place, rest mass energy is converted into internal energy 
and vice versa. In such a situation baryon number and 
lepton number are conserved, but not the total rest mass 
of the particles, and only the total ("relativistic") energy 
density, p + e, is a well defined quantity, but not the in- 
dividual energies. Case B therefore becomes ambiguous 
(see Sect. 12. 211 while Case A does not suffer from such 
problems. We therefore are tempted to recommend using 
Case A as an effective relativistic potential in Newtonian 
hydrodynamics codes. The extra factor T < 1 in the in- 
tegrand of the TOV mass, Eq. J3J, reduces the mass in- 



tegral below the mass equivalent of the total (rest mass 
plus internal) energy. Its introduction may be justified by 
heuristic arguments and consistency considerations when 
the relativistic concept of the gravitating mass is applied 
in a Newtonian description of the dynamics, which does 
not distinguish between proper volumes and coordinate 
volumes (see Sect. I2.3|) . 
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Appendix A: Construction of an "effective" 
one- para meter equation of state 

The relativistic TOV solution for a spherical, self- 
gravitating matter distribution in equilibrium can be 
found by solving the following two coupled ordinary dif- 
ferential equations for the pressure P and the mass m: 



dP 
dr 
dm 
dr 



p(l + e) + P , , N / 2to\ , k , 
(m + 4nr*P) (l-—), (AT) 

4?rp(l + e)r 2 . (A.2) 



with appropriate boundary conditions at the inner and 
outer boundary at r = and r = R, respectively. Note 
that for a TOV solution the mass m(R) is equivalent to the 
total gravitational mass M g of the matter configuration. 

The above system is closed by choosing a two- 
parameter EoS P(p, e) and assuming a distribution profile 
for the specific internal energy e, which effectively trans- 
forms the EoS into a one-parameter relation, P = P(p) [or 
alternatively p = p(P)], from which the specific internal 
energy can be determined as e = e(P, p). A typical exam- 
ple for this is to demand a polytropic relation for P and 
to specify the internal energy e by the ideal gas EoS: 



P = Kp~< 



P = pe(j-1) 



P= 77 



P 
K 
P 



l/ 7 



(A.3) 



P(7-1) 

For a two-parameter EoS of the form of Eqs. Q12I I13|) . 
where a priori no profile for the specific internal energy 
is known, the following procedure can be applied to con- 
struct an "effective" one-parameter EoS. For matter in 
equilibrium, e.g., at late times long after the core bounce 
phase, both the radial profiles of the pressure P(r) and 
the density p(r), as well as the specific internal energy 
e(r) can be assumed to be stationary. If the pressure pro- 
file is monotonic, then for each value of P there exists 
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a unique location r. By inverting the pressure profile, 
P(r) — ► r(P), we can thus relate a specific value p(r) 
and e(r) to each value of P. This way a one-parameter 
EoS p = p(r(P)) [and e = e(r(P))] can be constructed. 

The excellent matching of the density profiles in the 
central part of a collapsed core in equilibrium obtained 
from a dynamic evolution with profiles from a TOV solu- 
tion using equal central density and the EoS transformed 
as described above (see Figs. |21 to |SJ demonstrates the 
applicability of this method. 

Appendix B: Radial coordinate transformation 



To compare a solution of the TOV equations IjA.ll IA.2|) , 
which are formulated in standard Schwarzschild-like coor- 
dinates, to results from our evolution code using relativis- 
tic gravity, which is based on radial isotropic coordinates, 
we must perform a coordinate transformation of the radial 
coordinate. 

In standard Schwarzschild-like coordinates, the line el- 
ement ds 2 t reads 



ds s 2 t = -Adt 2 st 



Bdri 



(B.l) 



while the line element dsf so in radial isotropic coordinates 
is given by 



d 4o = dt L 



4 (dr? + r? 80 dfif sl 



(B.2) 



In vacuum outside the matter distribution, the metric 
components become 



A=l 
and 

1 



2M g 
r s t 



Ms 

2r iso 



B= 1- 



2AL 

& 



= A- 



1 + 



2r iso 



= 1 



Ms 



(B.3) 



(B.4) 



respectively. 

As with both metrics one can describe a spherical mat- 
ter distribution in equilibrium, equating the line elements, 
ds 2 t — ds 2 so , yields conditions for each of the coordinates 
and metric components in such a case. 

Setting the two coordinate times equal, we obtain a 
condition for A and a: 



b st ^ ^"iso' 

dt st = dtt s 



VA: 



(B.5) 



A similar condition can be obtained for the angle element: 

— ► r st = 2 r iso . (B.6) 



" 2 ; d^ 

dfLt — dfl\ 



2 

iso 7 



From the radial elements, we get a relation for the radial 
differentials: 



VB dr st = (j) 2 dr is 



Knowing the solution in the isotropic radial coordi- 
nate, it is straightforward to derive a simple expression 
for r st (r iso ) from Eq. (|B.6(I . 

On the other hand, obtaining the inverse relation 
nso( r st) is more complicated. For this, we insert Eq. (|B.6|) 
into Eq. (|B.7|I and arrive at 



(B.8 



dr st cj> 2 r st 

This differential equation can be solved by numerical 
integration using, e.g., a Runge-Kutta integration scheme 
with the following boundary conditions: 



7"0iso =r 0s t = 0, 
-Riso = — r~ I 1 — 



& 



2Mg 
-Rst 



(B.9) 
(B.10) 



An alternative method is to rewrite Eq. (|B.8p in terms 
of In r isn , 



(B.ll) 



(B.7) 



d In r iso _ 

dr st r st 
which can be directly integrated numerically. 
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